﻿ Metode de simulare a sistemelor analogice cu sisteme discrete Problema 1 Fie sistemul analogic liniar şi invariant în timp descris de funcţia de transfer: 2s+ Hs= () 2c 13ss++ ) ()( a) Să se determine răspunsul său la impuls, ( ; ) cht b) Să se determine răspunsul la impuls al sistemului numeric echivalent pe baza hn; metodei invarianţei răspunsului la impuls [ ] d Hz[hn ; c) Să se determine (), transformata z a lui d]d d) Să se obţină ecuaţia cu diferenţe finite pentru sistemul discret echivalent sistemului analogic considerat Rezolvare Punctul a) ABC + () 2cHs=+ 131sss++ + () 221 1s+ = =⇒= 1Bs HsB=+ ()() 11css=− − + 32 2s= 21s+ = =− ()() 23cCs Hs=+ 33ss=− 41s=− + () −− 1 2232ddsss++ ⎡⎤ ⎡⎤ = = = ()() 21cAsHs=+ ⎢⎥ ⎣⎦ 11ss=− 1s =− 343dsds ss=− + ⎣⎦ + () 11 1 1 11 +⋅ −⋅ () 2cHs=⋅ 412 431sss++ + () 77 311tt−− ete tσσ− ()() 44 11d⎛⎞ ⎜⎟ 2−= 11ds ss+ ⎝⎠ + () () ststdX s∞∞ −− L xt X sxt te dt txte txt s↔⇒ = − =− =− () ( )()( )() ()() {} ∫∫ −∞ ds−∞ 11 tt−− −=− ⇒−↔− ()()() {} 22Lte t ste tσσ 11ss++ () () 311 1ttt−−− =+ −σ ()()()() cht e t te t e tσσ 42 4 1 311nTnTnTT−−− Punctul b) [] == + −σ () [][][] dchn hnT e n ne n e nσσ 42 4 Punctul c) 11 1/4n nTT−− =↔ [][] () 1Ten e nσσ −− 44 1ez− 3311 1/4nnTT−− en e nσσ−=− ↔− [][] () 31T−− 44 1ez− 1T−− nTT Tez nTT−− ne n n e nσσ=↔ [][] () 12T−− ) 22 2(1ez− 1T−− 1/ 41/ 4Te z Hz=+ ⋅ − () 31T 112dTT−− −− − − 12(1)1ez ez e z−− − Punctul d) 211TTe− − −− − −−− 131 131 1TT T T−− 11 1 1ez e z z e z ez−− + −−− ()( )()() 424 () 2dHz== − − 131TT−− 11ez e z−− ()( ) − 44TT T−− 2T− − − − − −− − 31 2 1 21 21111 1TTTeTeTe−− eez z z z ez ez−+ + + − −+ − () 444 2 2 42 4 == − 44TTT−− ⎞ ⎛⎞⎛ −−−− 22131111TTTTeTeTe−− zezeee−− + −−+ ⎟ ⎜⎟⎜ 424 244 2 ⎠ ⎝⎠⎝ = 12231 4253TTT TT−− −− −− −− −− −+ − 122ez e z e z e z e z−+ 4423TTT TTTT−−− −−−− −− +−− 21222eTee eTeee−− zz+ 44 Hz= () 53 2 4 23 1dTTT TT−− − − − − − − 22ez e e z e e z−+ 1 + −+ + ()() 4422 31TTT TTT−−−−−−−− 22eTeezTeeez−− + +− ()() Hz= () 53 2 2 22 1dTTT TT−− − − − − − − 4+ 4412 42ez e e z e e z−+ + −+ () () 21−− 21bz bz+ = () − 321dHz−− + 0a 321az az az++ 22 22TT TTT−− −− − −=+−= ()() 0 2121; 21;0be e Te be T e b=− 522 2TTT TT−− − −− = + =− + = () () 32 1 04; 4 12 ; 4 2 ; 4aeae eae ea=− Ecuaţia cu diferenţe finite este : 2 522 2TTT TT−−− −− 43412 242 14eyn e eyn e eyn yn−−++ −−+ −+= [][][][] ()() 2222TT TTT−− −−− 21 1 2 1 2eT exn ee Te xn+− − + − − − [][] ()() Problema 2 Fie () funcţia de sistem a unui sistem analogic liniar şi invariant în cHs 1 Hs= timp, () 2c 12ss++ ()( ) se determine răspunsul la impuls al sistemului analogic, () ; a) Să cht b) Se notează cu () răspunsul indicial al sistemului analogic Se notează cu cst hn răsn ră spunsul la impuls şi cu [spunsul indicial al unui sistem discret []] dd liniar şi invariant îin timp În ipoteza că sistemul discret s-a obţinut prin echivalarea sistemului analogic considerat folosind metoda de echivalare bazată , este adevarată relaţia : pe invarianţa răspunsului la impuls, adică [ () ] cdhn hnT= n ? () [] ∑ cdsn hkT= k=− ∞ sistemul discret echivalent s-a obţinut prin metoda invarianţei răspunsului c) Dacă , este adevarată relaţia :[ ? indicial, adică [ ()() ]] cdsn snT=cdhn hnT= ei răspunsului indicial, să se determine d) In ipoteza de la punctul c), a invarianţ Hz () pentru sistemul discret echivalent d Rezolvare Punctul a) 1ABC ++ () 22cHs== 1212 2ssss s++ ++ + ()( ) () 1 = = ()() 211cAs Hs=+ =− 11ss=− 2s+ () 21 = =− 21Bs Hs=+ ()() 22css=− + 1s=− 211dd⎛⎞ ⎡⎤ − = =− = ()() ⎜⎟ 221cCsHs=+ ⎣⎦ =− 2 22sss=− 11dsds ss=− + ⎝⎠ + () 22ttt−−− σ=− − => () ()()() cht e t te t e tσσ Punctul b) nn∞ =∗ = −= =kT () [] [] [] [][ ] [] ∑∑∑ cddd dsnnhnhknkhkhσσ kkk=−∞=−∞ =−∞ ∞ Deci : [] Relaţia propusă este adevarată în cazul general [ ] ∑ cdsn hkT= ∞ k=− Punctul c) Fie ( transformata Laplace a raspunsului indicial al sistemului analogic ) cSs 3 () 1cHsAB C D =++ + () 2cSs== + + 2 12 12ssss ss ss++ ()( ) + () Se calculează răspunsul indicial al sistemului analogic : 11 AsSs=== () 00css== 12ss++ 4 ()( ) 1 = =− ()() 211cBs Ss=+ =− 11ss=− 2ss+ () 211 2Cs Ss=+ = = ()() 22css=− − 12ss= + () 21213dds+− ⎡⎤ 23DsSs⎡⎤ /4 =+ = =− =−= ()() 22c⎢⎥ ⎣⎦ =− 222sss=− 4dsds s sss=− + + ⎣⎦ () 22113ttt−− − =>() σ σ=− + +t () ()()() cst teteteσσ 424 Raspunsul la impuls al sistemului echivalent este: 22nTnTnT=−− == − −σ () [][][][] dchn hnt e n nTe n e nσσ sn snT=⇒ () [] cd 22113nTnTnT−−− sn n e n nTe n e nσσ=− +⋅ + σσ [][] [][][] d 424 21T−− 11 13 1Tez⋅ Sz=⋅ − +⋅+⋅ () 11 22dTT1−−−− − 21424111Tzezez−− −−⋅−⋅ 1ez−⋅ () 22 − − − − − 121 121TT T−− 11 411ez e z z e z−⋅ −⋅ −− −⋅ ()( )()() Sz= () 2d −−− 1121TT−− 41 1 1zezez−− ⋅−⋅ ()( )( ) − −− − −− −− 1T 2111 112TT T−− 211 311 1Tezzez zezez⋅− −⋅ − −⋅ −⋅ ()( )()( )() ++ 22 −−− −−−−− 1 1121112TT TT−− 41 1 141 1 1zezez zezez−− ⋅−⋅ −−⋅−⋅ ()( )( )()( )() 1− 1Hz z Sz=− = ()() () dd 12142 12142TTT TT−− −−−− − −−−− 112 4112ez ezez z ezez−⋅ −⋅+⋅ −− −⋅+⋅ ()( )()() + =+ 22 121121TT TT−− − −−− − − 41 141 1ez e zez e z−⋅ − ⋅−⋅ − ⋅ ()( )()( ) 212 1 1 1 2 21TT TT−−− −− −−−−− −− T 2131 1Te z z e z z e z e z e z−−⋅ −−⋅+⋅−⋅ ()( )()() ++ 22 121TT−− − −− 121TT−− − 1zez⋅−⋅ 41 141ez e ze−⋅ − ⋅− ()( ) ()( ) 4 2142 1 3253TTT TT−−−−−− −−−− 122ezezez ezez⋅⋅+⋅⋅⋅+ − Hz= () 2d 121TT−− − − 41 1ez e z−⋅ − ⋅ ()( ) 42 21 143 22TT TT−− −− −−− −− 41 22ez ez zez ez−+ − − − () + 2 121TT−− − − 41 1ez e z−⋅ − ⋅ ()( ) 21 22 3TT T−− −−−−− 2Te z e z z e z−−+ () + 2 121TT−− − − 41 1ez e z−⋅ − ⋅ ()( ) 11TT−− 221223233TTTT−−−−−−−− − − − 31ze z e z−+ ez ez ez ez−−+− () + 2 121TT−− − − − ⋅ 41 1ez e z−⋅ ()( ) 323TTT−− − − +− 423ee eTz−+ () = 2 121TT−− − − − ⋅ 41 1ez e z−⋅ ()( ) 24 3 2TTTT−− − − − −++ +− + 352 52 3z e Te Te e⎡⎤ () () ⎣⎦ + 2 121TT−− − − − ⋅ 41 1ez e z−⋅ ()( ) 312TT T−− − − −+++ 524 52 1Te z e Te++ ()() () + 2 1T−− 21T−− − 41 1ez e−⋅ z⋅ () () Deci: 32343TT TTTT−− − − 2−−−− 423352 3ee e TzeTeez⎡⎤ ++− −++ + () () ⎣⎦ Hz=+ () 22d − −−− − − 121121TT TT−− 41 141 1ez e z ez e z−⋅ − ⋅ −⋅ − ⋅ ()( ) ()() − 21TT−− 452 1eTez−+ + + () () + 2 − − 121TT−− 41 1ez e z−⋅ − ⋅ ()( ) Adică: 5 32332TT TTT T T−− − − 2−−−−− 423352 52 3ee e Tze e Te Te z⎡⎤ ++− −++ +− + () () () ⎣⎦ Hz=+ () 22d − −−−− −− 121121TT TT− 41 141 1ez e zez e z−⋅ − ⋅−⋅ − ⋅ ()( ) ()( ) α 21TT−− − 452 1eTez⎡⎤ −++ + () ⎣⎦ + 2 − − 121TT−− 41 1ez e z−⋅ − ⋅ ()( ) AB C + Hz=+ () 1221dTT−−− − 2111Teze z−− −⋅− ⋅ 1ez−⋅ () 1T−− =− ⋅ () () 11TdAezHz− ze= − − 323TTT−− +− 423ee eTz−+ () = 2 21T−− 41ez−⋅ () 1T− ze= − − −− 32 21TT T T−− 352 52 32e e Te Te z z⎡⎤ −++ +− + + () () ⎣⎦ + 2 21T−− 41ez−⋅ () 1T− ze= 87 62TT TT−− −−T 42334ee Tee−− +− − +e () A= 2 T− 41e− () nnn 2TT−−T− 1hn Ae n Bn e n Ce nσσ=++ +σ () [][][][] ()() () d 22nTnTnT−−− σ=− − () [] [] [] chnT e n nTe n e nσσ Ahn hnT≠⇒ ≠1 Deoarece: [ ia propusă nu este adevarată () ] cddeci relaţ Problema 3 Se notează () răspunsul la impuls al unui sistem liniar, invariant în timp cht hn cel al unui sistem liniar, invariant în timp discret ce echivalează continuu şi cu [ ] d sistemul analogic prin metoda invarianţei răspunsului la impuls t− = () să se determine functia de transfer a filtrului analogic şi să a) Dacă () cht e tσ se deseneze dependenta de frecvenţă a modulului său ; se determine funcţia de transfer a filtrului numeric ; să se schiţeze b) Să caracteristica de modul a răspunsului în frecvenţă al sistemului discret c) Sa se repete punctul b) astfel încât caracteristica filtrului numeric sa fie aproximativ identica cu cea a filtrului analogic Rezolvare Punctul a) 6 111 t− Hω ⇒=⇒= ==L ()() () () {} cHωc cHs e tσ 2 1jω 1s+ + 1ω + Deci sistemul este de banda nelimitata => apare alierea la frecvente mari Pentru T=0 1: 11ππ ⎛⎞⎛⎞ 0 1; 0, 063; 0, 031HHH== ≅ == () ⎜⎟ 22ccc⎜⎟ 2TT⎝⎠ ⎠ ππ 1(5)1(10)⎝ ++ Punctul b) 11 nT− hn hnT e n H zHσ == ⇒= ⇒Ω= ()()() [][] 1dddjTTc−Ω −− − 11ezee− −⋅ 11 HΩ== () 2dTT−− −⋅Ω+ ⋅Ω 221cos sinTTeje−− ⋅Ω+ ⋅ Ω 1cos sinee− () 1 = 2TT−− Ω 12cosee+− 111 1π1 ⎛⎞ HHHπ0;; =≅ = = = ()() ⎜⎟ d 2112ddTTTTT−− () 11− ee−− ⎝⎠ 1e−+ + Se ia perioada de eşantionare, T= 0 1, astfel incat frecventa f =1/T, sa fie de cel puţin 10 ori mai mare decât frecventa naturala a sistemului, inversul constantei sale de timp 0,20,1−− 0 10, 50; 2 1 1 0, 74; 1 1 0, 52HH eHeππ==+≅=+= ()()() () dd d 7 Punctul c) Pentru ca spectrul sa fie aproximativ identic pentru Ω ∈ (-π,π) respectiv ω ∈ (-π/T,π/T) (se stie Ω = ωT) : nTT− h n Th nT Te n H zσ== ⇒= ()() [][] 1ddTc−− 1ez− TT HHΩ=⇒ Ω = ()() − 2ddjTTT−Ω 1ee−−−⋅ ees+ 12co−Ω Pentru T=0 1 : 0,20,1−− 0 1, 05; 2 1 1 0, 074; 1 1 0, 052HH e Heππ==+≅=+= ()()() () dd d 8 Problema 4 Se consideră sistemul analogic liniar şi invariant în timp descris de ecuaţia: dy () 2yt xt+= () dt a) Să se determine răspunsul la impuls al unui sistem discret ce simulează sistemul dat prin metoda invarianţei răspunsului la impuls, precum şi funcţia sa de sistem ; se determine aceleaşi funcţii ca la punctul a) pentru un sistem discret b) Să echivalent, obţinut prin metoda aproximării numerice a ecuaţiei diferenţiale ; se repete punctul a) pentru un sistem discret obţinut prin metoda transformării c) Să biliniare aţi formele de implementare pentru sistemele de la punctele a), b) şi c) ; d) Schiţ sinxt=t ş i o perioada de esantionare 0,1T= să se determine e) Pentru () amplitudinile semnalelor de ieşire ale sistemului analogic dat şi ale celor trei sisteme numerice Rezolvare Punctul a) () 1Ysdy 2 yt xt sY s Y s X s H s+= 2 + = ⇒ == () () () () () () a⇒ 2dtX s s+ () 221tnT−− ht e t hn e n H zσσ =⇒= ⇒= () ()() [] [] 21addt−− 1ez− Punctul b) Aproximarea numerica a ecuaţiei diferenţiale se face prin aproximarea derivatei: 1− 11yn yndy tz−− () [][] − s=⇒= dtTT tnT= 1− 112yn ynTz−− [][] +− 2ynxnYz Xz⇒+=↔ = () () [] [] dd TT T Hz⇒= () 1d− 12Tz+− Sau se înlocuieşte direct în transformata Laplace: T 1T + 112THz Hs− () ==== () − 111zads− − = − +− 11112TzTz− ⋅ 12z− + 12TT+ n 1T⎛⎞ = hnnσ [][] d⎜⎟ 12 12TT++ ⎝⎠ 9 Punctul c) Aproximarea biliniara: 1− 1Tz+ () 1 1− === () () − 121zadsHz Hs− = 11−− − 12121 1TzzT z− −+ + () () 1z+ + 12− 1Tz+ 1− 1− 1Tz+ () Tz+ 1 ==⋅ 11−− 21 21 2 1 1TTz TTz+− − +−− ()()() ⎡+ +⎤ ⎛⎞ 111TT− 11Tz−− ++ () ⎜⎟ ⎢⎥ 11TTT−− ⎝⎠ ⎣⎦ =⋅ 1− 211 1TT Tz⎡⎤ −+−− () () ⎣⎦ + ⎛⎞ 11T− 1T+ 1zT−− () + 1⎜⎟ 1TTT− ⎝⎠ 1T− =⋅ +⋅ 1− + 11TTz+−− − ⎛⎞ ()() 11212TT− 1Tz−− () ⎜⎟ 1T− ⎝⎠ 21TT− =+⋅⋅ 1− 21 2 1 1 1TTTTz−− +−− ()() 1TT =− + ⋅ 1− 21 1 1 1TTTTz−− +−− ()() 1TT =− +⋅ − −−+ 1121 1 1TTTT− ()()() 1z− 1T+ n 1TTT− ⎛⎞ hnnnδσ =−+ ⋅ [][][] 2d⎜⎟ 211TT−+ () 1⎝⎠ T− () Punctul d) x y Pentru sistemul analogic ∫ 2yyx+= ∫∫ -2 ∫ 10 x[n] y[n] Pentru sistemul numeric de la a) 2T− baae===−1; 1; 001 -2TD 2T− 1yn e yn xn−−= [ ] [] e[] -1 x[n] (1+2T)y[n] Pentru sistemul numeric de la b) T ; 1 2 ; 1bTa Ta==+ =− 00 1 121Tyn yn Txn+−−= [] [ ] [] D () 1 y[n-1] Pentru sistemul numeric de la c) 1 bb==1; 1; 01 (y[n]x[n] 21T+ ) 21 ; 21aTaT=+ =− () () 01 +− −= 21 21 1Tyn Tyn+ () [][] D D () − 1xn xn=+ [] [ ] -2(1-T) Punctul e) sin ; 1 sin arg 1xt t yt H t H== + () () ()() {} () 1111 1 1Hs HHHω=⇒= ⇒=⇒= ()()()() 22 25sjjω++ + 11⎛⎞1 arg 1 arg 2 arctg sin arctgHjytt=− + =− ⇒ = − (){}() {} ⎜⎟ 225⎝⎠ 1 yA⇒= 5 Pentru sistemele numerice: 11nn⎛⎞⎧⎫ ⎛⎞ ⎛⎞⎛⎞ sin 0,1 sinsin argxn nyn HH==⇒= + () [][] ⎬ ⎜⎟ ⎜⎟ ⎜⎟⎜⎟ dd⎨ 1010 10 10⎝⎠ ⎝⎠⎝⎠ ⎩⎭ ⎝⎠ Pentru sistemul numeric de la a) 111 () Hz=H⇒Ω= = () −− 21dT−−222djTTT−Ω − Ω+ ⋅ Ω 1ez−⋅ 1cos sin1ejeee− −⋅ 11 11 H⇒Ω== () 24dTT−−2 24212coTTee−− −Ω s+ + Ω 1cos sinee−Ω () 11 1⎛⎞ H=≅ 2d⎜⎟ 0,2101− ⎝⎠ 0,4 0,2−− 1e− ee+− 12cos () 10 Pentru sistemul numeric de la b) TTT HzH=⇒ Ω= = ()() − 1ddj−Ω 12cos sin1212TjTzTe+Ω− Ω +⋅ +⋅ TT H⇒Ω == () 222d 14 4cos12cos 2sinTTTT++ Ω +Ω+ Ω ()() 10,10,1⎛⎞ 0, 083H⇒≅== d⎜⎟ 101, 21, 44⎝⎠ Pentru sistemul numeric de la c) Ω − 1j− 1111ze++ HzH=⋅⇒ Ω=⋅ ()() − 1ddj−Ω 22111, 1 0, 9TTze+− − − () 22 1cos sin1122cos+Ω + Ω () +Ω =⋅=⋅ H⇒Ω () 22d 22 1, 21 0,81 1,98cos+ 1,1 0,9cos 0,9sin−Ω −Ω+ Ω ()() valoarea 0,1 si se calculează Se dă lui Ω Problema 5 Fie un sistem în timp continuu cu răspunsul la impuls ()ht şi funcţia de c transfer () Ecuaţia diferenţială care descrie acest sistem este: cHs kkNN dy dx ab= kkkk∑∑ 00kkdt dt== Se aproximează acest sistem cu unul în timp discret folosind aproximarea: 1xnxdx+−n [][] ≅ tnTdtT= 0 ) ( Se definesc diferenţele de ordinul 0 : [] xnxn∇= [] {} 11xnx+− n [][] ) ( de ordinul 1 :[] xn∇= {} T 1kk− ()() 1 xnx∇=∇∇n şi de ordinul k :[] [] {}{} {} 12 s Să se determine funcţia de transfer a sistemului discret a) Fie () Hs= 2c 1s+ () Hz echivalent () ; d tura dintre (Hs şi (Hz ? b) Care este legă )) c d Rezolvare 11xn xn+− [][] ) ( xn∇= [] {} T 11xn xn⎧⎫ [][] ) ( 2+− xn∇= ∇ [] {} ⎨⎬ T ⎩⎭ 111⎡⎤ )) (( 1xn xn=∇ +−∇ [][] {}{} ⎢⎥ T ⎣⎦ 2111xnxnxnxn⎡⎤ +− + +− [][][][] =− ⎢⎥ TT T ⎣⎦ 1 2 2 1xn xn xn⎡⎤ =+−++ [][][] 2⎣⎦ T Se observă că: k 1l kl +− 1xn Cxn l∇= () [][] {} kk∑ 0lT= Punctul a) bss 1 Hs== () 22c + + 21ss aasas++ 01 2 = = = = 0; 1; 1; 2; 1bbaaa⇒= 0101 2 Ecuatia diferenţială corespunzătoare este : 2 dy t d y t dx t () ()() 2yt++ = () 2 dtdtdt Ecuaţia cu diferenţe finite corespunzătoare este : 121 ))) ((( 2yn yn yn xn+∇ +∇ =∇ [][][][] {}{}{} Adică : 111ynynxnx+−+−n ][][][] [ +++−++= 2221ynyn yn yn⎡⎤ [][][][] 2⎣⎦ TTT Sau : 2 21 221 1T n T n n n n n Tn Tnyyyy yyx++−++−++=+−x [][][][][][][][] () 13 Luând în ambii membrii ai ultimei relaţii transformata Z, se obţine: 22 212TzTz zzzzz zTzzTzYYYYYX+− + − += −X () ( ) () () () () () () Sau : 22 22 21Yz T Tz T z z Xz Tz T⎡⎤+−+−+= − ()()( ) ⎣⎦ Tz T− Hz⇒= () 22d +−+ 21 21zTzTT+− () Punctul b) 1TzTz T− () − Hz== () 22d +− 11zT zT+− () () 1z− 11Tz− () T ⋅ T== 22 T ⎛⎞ ⎛⎞ 211zz−− 11T++ ⎜⎟ ⎜⎟ TT⎝⎠ ⎝⎠ 1z− T Deci : () Hz= 2d 1z− ⎛⎞ 1+ ⎜⎟ T⎝⎠ Se constată că : () = () 1zcHz Hs− ds= T 14